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Abstract — We study throughput-optimum localized link 
scheduling in wireless networks. The majority of results on 
link scheduling assume binary interference models that simplify 
interference constraints in actual wireless communication. While 
the physical interference model reflects the physical reality 
more precisely, the problem becomes notoriously harder under 
the physical interference model. There have been just a few 
existing results on link scheduling under the physical interference 
model, and even fewer on more practical distributed or localized 
scheduling. In this paper, we tackle the challenges of localized 
link scheduling posed by the complex physical interference 
constraints. By cooperating the partition and shifting strategies 
into the pick-and-compare scheme, we present a class of local- 
ized scheduling algorithms with provable throughput guarantee 
subject to physical interference constraints. The algorithm in 
the linear power setting is the first localized algorithm that 
achieves at least a constant fraction of the optimal capacity 
region subject to physical interference constraints. The algorithm 
in the uniform power setting is the first localized algorithm 
with a logarithmic approximation ratio to the optimal solution. 
Our extensive simulation results demonstrate correctness and 
performance efficiency of our algorithms. 

Index Terms — localized link scheduling, physical interference 
model, maximum weighted independent set of links (MWISL), 
capacity region. 

I. Introduction 

AS a fundamental problem in wireless networks, link 
scheduling is crucial for improving networking perfor- 
mances through maximizing throughput and fairness. It has 
recently regained much interest from networking research 
community because of wide deployment of multihop wire- 
less networks, e.g., wireless sensor networks for monitoring 
physical or environment [1| [2| through collection of sensing 
data O. Generally, link scheduling involves determination of 
which links should transmit at what times, what modulation 
and coding schemes to use, and at what transmission power 
levels should communication take place |4|. In addition to 
its great significance in wireless networks, developing an 
efficient scheduling algorithm is extremely difficult due to 
the intrinsically complex interference among simultaneously 
transmitting links in the network. 
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The link scheduling problem has been studied with different 
optimization objectives, e.g., throughput-optimum scheduling, 
minimum length scheduling. Our study mainly focuses on 
maximizing throughput in multihop wireless networks. Nev- 
ertheless, a common primary issue here is to find a set of 
conflict-free links with regard to the interference models. 

Interference models matter significantly in scheduling al- 
gorithm design and analysis. It is well known that the time 
complexity of a throughput-optimum scheduling that tries 
to find a maximum weighted independent set of links, is 
polynomial under the 1-hop interference model in switched 
wired networks, whereas it becomes NP-hard under more 
general interference models in wireless networks fT]. Here 
the weight associated with every link is its corresponding 
queue length. Motivated by providing efficient solutions for 
the NP-hard problem, numerous low complexity scheduling 
algorithms with suboptimal throughput guarantee Q, H, ISl, 
0, H, 0, d, ini, 113, lHai, have been proposed in 
literature. 

Despite of the numerous significant results gained for the 
problem, most of them assume simple binary interference 
models, e.g., hop-based, range-based, and protocol interference 
models [14J. Under this category of interference model, a 
set of links are conflict-free if they are pairwise conflict-free. 
Conflict of transmissions on two distinct links is predetermined 
independently of the concurrent transmissions of other links. 
Thus, interference relationships based on these models can be 
represented by conflict graphs, and we can leverage classic 
graph-theoretical tools for solutions. However, in actual wire- 
less communication, interference constraints among concur- 
rent transmissions are not local and pairwise, but global and 
additive. Conflict of distinct transmissions is determined by 
the cumulative interference from all concurrent transmissions, 
which is often depicted by the physical interference model, 
e.g., the Signal-to-Interference-plus-Noise Ratio (SINR) inter- 
ference model. But due to the simplification of binary inter- 
ference models on physical reality, the corresponding results 
no longer hold under the more realistic physical interference 
model. Recent results lITSl . |[T6ll have shown that throughput 
differs significantly under the two kinds of models. Moreover, 
the global and additive nature of the physical interference 
model drives previous traditional techniques based on conflict 
graphs inapplicable or trivial. Consequently, designing and an- 
alyzing scheduling algorithms under the physical interference 
model becomes especially challenging. 
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Some recent research results [HT], IH, [HI], [HI], [[20], 
[El, [US, [113, [124 1 have addressed a few challenges related 
to the link scheduling problem under the physical interference 
model. To the best of our knowledge, however, all of these, 
with throughput maximization or other optimization objectives 
such as a minimum length schedule [14], just focus on central- 
ized implementation. Distributed or even localized scheduling 
under the physical interference model is more demanding out 
of practical relevance. 

Though [|25]| , [1261 , [1771 consider distributed implementation 
of centralized algorithms, they fail to provide an effective 
localized scheduling algorithm with satisfactory theoretical 
guarantee [|25l [|26l and require global propagation of mes- 
sages. A scheduling algorithm without theoretical guarantee 
may cause arbitrarily bad throughput performance, and global 
propagation of messages throughout network is inefficient in 
terms of time complexity [28 ], especially for large-scale net- 
works. This motivates us to develop localized link scheduling 
algorithms with provable throughput performance. Here by a 
localized scheduling algorithm we mean that each node only 
needs information within constant distance to make scheduling 
decisions, while a distributed algorithm may inexplicitly need 
information far away. Since just local information is available 
for each node to collaborate on schedulings with globally 
coupled interference constraint, it poses significant challenges 
to designing efficient localized scheduling algorithms with 
theoretical throughput guarantee. 

In this paper we tackle these challenges of practical lo- 
calized scheduling for throughput maximization under the 
realistic physical interference model with the commonly- 
used linear and uniform power assignment. For the primary 
challenge of decoupling the global interference constraints, 
on observing that distance dominates the interference, we 
partition the links into disjoint local link sets with a certain 
distance away from each other. By this way, we bound the 
interference from all the other local link sets with a constant 
so that independent scheduling inside each local link set is 
possible. Then we shift the partitions to ensure every link 
get scheduled. Using this we implement localized schedulings. 
Furthermore, we novelly combine these two strategies into the 
pick-and-compare scheme by which we can prove that our 
algorithms obtain a theoretical capacity guarantee. 

We summarize our main contributions as follows. 

1) We propose a localized link scheduling algorithm under 
the linear power setting with provable throughput guar- 
antee. We successfully decouple the global interference 
constraints for the implementation of localized schedul- 
ing, by proving that there exists a constant distance 
between any two disjoint local link sets to ensure inde- 
pendent link scheduling inside a local link set under the 
linear power setting. We then prove that our algorithm 
can achieve a constant fraction of capacity region. 

2) We further revise the algorithm to work under the 
uniform power setting, and give a detailed analysis on 
the achievable capacity region. While under the uniform 
power setting, we prove that if the size of the scheduling 
set for local link sets is upper bounded by a constant, the 
solution to decouple the global interference constraints 



is also applicable. We prove that our revised algorithm 
can achieve 0(log \ V\) fraction of capacity region where 
\V\ is the number of nodes. 
3) We conduct extensive simulations to evaluate our al- 
gorithms. The simulation results demonstrate the cor- 
rectness of our algorithms and corroborate our theo- 
retical analysis. The comparisons with two centralized 
algorithms in terms of average throughput performance 
further show performance efficiency of our algorithms. 

The remainder of the paper is organized as follows. In 
Section II, we define the exact system models for our problem. 
In Section III, we describe the basic ideas of our solutions, and 
the proposed localized scheduling algorithm with theoretical 
performance analysis. In Section IV, we revise our solutions 
for the uniform power setting, and propose a localized algo- 
rithm with a logarithmic approximation ratio. We evaluate our 
solutions for both power settings in Section V. In Section VI, 
we review related works on link scheduling. In Section VII 
we discuss some related issues on proposed algorithms. We 
conclude this paper in Section VIII. 

II. Models and Assumptions 
A. Network Communication Model 

We model a wireless network by a two-tuple (V^E), where 
V denotes the set of nodes and E denotes the set of links. 
Each directed link / = (u^v) G E represents a communication 
request from a sender u to 3. receiver v. Let ||/|| or \\uv\\ 
denote the length of link /. We assume each node knows its 
own location, which is necessary in the partition and shifting 
process. 



B. Interference model 

Under the physical interference model, a feasible schedule 
is defined as an independent set of links (ISL), each satisfying 
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where ^ denotes the ambient noise, a denotes certain thresh- 
old, and Tu denotes the set of simultaneous transmitters with 
u. It assumes path gain r] ■ \\uv\\~'^ < 1, where the constant 
K, > 2 is path-loss exponent, and r] is the reference loss factor. 
We consider the following two transmission power settings. 

1) Linear power setting: a sender u transmits to a receiver 
V always at the power = c ■ \\uv\f < P where c 
and /3 are both constant satisfying c > 0^0 < /3 < hz, 
and P is the maximum transmission power. 

2) Uniform power setting: all links always transmit at the 
same power Pu = P. 

We use Pi and Pu alternatively to denote the transmitting 
power of link / = {u^ v). 

We also assume all links having a length less than the 
maximum transmission radius The distance between 

u and V satisfies r < \\uv\\ < R, where r and R respectively 
denotes the shortest link length and the longest link length. 
We suppose that r and R are known by each node. 
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C. Traffic models and scheduling 

The maximum throughput scheduHng is often studied in the 
following models. It assumes time-slotted wireless systems, 
and single-hop flows with stationary stochastic packet arrival 
process at an average arrival rate A^. The vector A{t) = 
{Ai{t)} denotes the number of packets arriving at each link 
in time slot t. Every packet arrival process Ai{t) is assumed 
to be i.i.d over time. We also assume all packet arrival process 
Ai{t) have bounded second moments and they are bound by 
Amax, i.e., Ai{t) < Amax,V/ G E. Let a vector {0,1}I^I 
denote a schedule S{t) at each time slot t, where Si{t) = 1 if 
link / is active in time slot t and Si (t) = otherwise. Packets 
departure transmitters of activated links at the end of time 
slots. Then, the queue length (it is also referred to as weight 
or backlog) vector Q{t) = {Qi{t)} evolves according to the 
following dynamics: 

Q{t + 1) = max{0, Q{t) - S{t)} + A{t + 1). 

The throughput performance of link scheduling algorithms 
is measured by a set of supportable arrival rate vectors, named 
capacity region or throughput capacity. That is, a scheduling 
policy is stable, if for any arrival rate vector in its capacity 
region llT2l . 

lim E[Q{t)] < oo. 

A throughput-optimal scheduling algorithm can achieve the 
optimal capacity region. It is well known that the policy 
of finding a maximum independent set of links to schedule 
regarding to the underlying interference models is throughput- 
optimal ||29l. Unfortunately, the problem of finding a MWISL 
itself is NP-hard generally [|3Ql . Thus we have to rely on 
approximation or heuristic methods to develop suboptimal 
scheduling algorithms running in polynomial time. A subopti- 
mal scheduling policy can just achieve a fraction of the optimal 
capacity region. The fractionally guaranteed capacity region is 
depicted by efficiency ratio 7 EI. 

A suboptimal scheduling policy with efficiency ratio 7 must 
find a 7-approximation scheduling at every time slot t to 
achieve 7 times of the optimal capacity region [ 31J . It remains 
difficult to achieve in a decentralized manner. The pick-and- 
compare approach proposed in [5 1 enables that we just need 
to find a 7-approximation scheduling with a constant positive 
probability. The basic pick-and-compare [5 | works as follows. 
At every time slot, it generates a feasible schedule that has a 
constant probability of achieving the optimal capacity region. 
If the weight of this new solution is greater than the current 
solution, it replaces the current one. Using this approach 
achieves the optimal capacity region. The proposition below 
further extends this approach to suboptimal cases. 

Proposition 1: ( flU]) Given any 7 G (0, 1], suppose that 
an algorithm has a probability at least S > of generating an 
independent set X{t) of links with weight at least 7 times the 
weight of the optimal. Then, capacity 7 • A can be achieved by 
switching links to the new independent set when its weight is 
larger than the previous one( otherwise, previous set of links 
will be kept for scheduling). The algorithm should generate 
the new scheduling S{t) from the old scheduling S{t — 1) 
and current queue length Q(t). 



In the rest of the paper we look for solutions to generate 
localized link schedulings of approximately optimal weight, 
with a certain constant probability at least. 

III. The Algorithm in the Linear Power Setting 

In this section, we focus on the design of localized schedul- 
ing algorithm for the linear power setting. We firstly demon- 
strate the basic ideas of the algorithm design before involved 
in the details of implementation. 

A. Basic idea 

The basis of our idea is to create a set of disjoint local 
link sets in which the scheduling can be done independently 
without violating the global interference constraint. The de- 
coupling of the global interference constraint is based on 
the fact that distance dominates the interference between two 
distinct links. That is, if a transmitting link is placed a certain 
distance away from all the other transmitting links, the total 
interference it receives may get bounded. We prove later that 
this assumption holds actually. Based on this, then we employ 
the partition strategies to divide the network graph into disjoint 
local areas such that each local area is separated away by 
a certain distance to enable independent local computation 
of schedulings inside every local area. Links lying outside 
local areas will keep silent to ensure separation of local areas. 
As links lying outside local areas cannot remain unscheduled 
all the time or it will induce network instability, we use the 
shifting strategy to change partitions at every time slot to 
make sure that every backlogged link will be scheduled. These 
locally computed scheduling link sets compose a new global 
schedule X{t) at every time slot t. 

In light of the pick-and-compare scheme, we choose a more 
weighted schedule, denoted as S{t), between a new generated 
schedule X(t) and the last- time schedule S(t — 1) using Q(t). 
Meanwhile, by Proposition 1, if we guarantee that 

f>{S{t)-Q{t)>^S%t)-Q(t))>5 (1) 

for some constant 7 > 0, 5 > 0, the queue length vector 
Q{t) will eventually converge to a stable state. Then we can 
get a constant approximation ratio scheduling policy for the 
optimal. 

B. Detailed description 

Now we describe the details of our method. 

We first give a brief description of the partition and shifting 
strategies ifTTI . as illustrated in Fig. L We partition the plane 
into cells with side length d = R, using horizontal lines x = i 
and vertical lines y = j for all integers i and j. A vertical strip 
with index i is {{x^y)\i < x < i + 1}. Similarly, we define the 
horizontal strip j. Let cell(i, j) denote the intersection area of a 
vertical strip i and a horizontal strip j. A super-subSquare(i, j) 
is the set of cells: {ce//(x, y)\x ^ [i * i^T + a^, (i + 1) * i^T + 
at),y ^ [j*i^ + 6t, + + and a sub-square(i,j) 

inside it is the set of cells: {ce//(x, y)\x e [i ^ K -\- at -\- M , {i -\- 
1) ^ K ^ at - M),y e [j * + 6^ + M, (j + 1) * + 
bt — M)}. The corresponding link set Yij (or Lij) consists 
of links with both ends inside super- sub Square(i,j) (or sub- 
square(z, j)). Let constant K = 2M + J, < at^h < K. 
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TABLE I: Summary of notations 



J 


side length of sub-square 


Lij 


link set of sub-square (i, j) 




link set of super-subSquare(i, j) 


K 


side length of super-subSquare 




new scheduling for L^j at t 


Sij(t) 


scheduling for Yij at t 


Zi 


local link set 




local optimal MWISL for L^j 


S*(t) 


global optimal MWISL at t 


R 


longest link length 


Shit) 


intersection of Lij and S*(t) 


S(t) 


global scheduling at time slot t 


d 


side length of cell 


\\uv\\ 


link length 


rs(l) 


relative interference / get from link set S 


asil) 


affectness I get from link set S 


Q{t) 


queue length vector 


W{S) 


weight of link set S 




interference link / suffered from link set S 




the maximum interference / can bear 


Imax 


the minimum of I in ax 



Super-subSquare Sub-square 



^3 



/ > 



Super-subSquare Sub-square 




(a) Partition (K, at, bt). Here the 
gray area is the local area the 
link set of which participate in 
computing the new schedule; the 
links in the white area keep 
silent. 



(b) Partition(ir, a(t+i), 6t). 
Change the partition through 
shifting to the right by one cell, 
ensuring that links in the white 
area in previous partitions have 
opportunity to be scheduled. 



Fig. 1: The partition and shifting process. The partition 
changes at different time slots by shifting vertically or hor- 
izontally by a step of one cell-length. 

at, bt are adjustable variables, referred to as the horizontal 
and vertical shifting respectively. M is a constant, relating 
to the distance for independent scheduling in local areas. We 
will formally define it later in Lemma 2. We define the above 
process as Partition{K^at^bt). By separately adjusting at, bt 
by a step of single cell, we can get K'^ different partitions for 
a plane totally. 

Note that herein nodes do not shift, it looks like that 
a virtually partitioned cover having the same size with the 
original topology shifts vertically or horizontally. Initially, the 
virtual cover and the original topology of the network coincide 
completely. When the virtual cover shifts, nodes located in the 
edge area will certainly get uncovered if the distance between 
the edge of topology and the nearest nodes is less than Kd. 
To tackle this, we just need to assume a larger virtual cover, 
e.g., the distance between the edge of topology and the nearest 
nodes is equal or greater than Kd. 

Then at each time slot nodes cooperate to compute a 
globally feasible scheduling as follows. 

Time slot t = 0: Every node first decides in which 
cell it resides by a partition using (ao,6o) = (0,0); then it 
participates in the process of computing a local scheduling 
Sij{0) for the sub- square it belongs to. Let the solution 
5'(0) of time slot be the union of the local solutions Sij{0) 
for all sub- squares. 

Time slot t > 1: Every node decides in which cell it resides 
by a partition starting from (at^bt). The shifting strategy for 
at and bt works as follows. We let at = t mod K; and bt = 
{bt + 1) mod K if at = 0, or it keeps unchanged. Each 



node then participates in computing the new local scheduling, 
denoted as ^ij{t), for its sub- square using the weight 
Q{t). Let Sij (t- 1) be the set of links from l)(the global 
solution at time slot t — 1) falling in the super-subSquare(z, j) 
instead of sub-square(i, j). If Sij{t — 1) • Q{t) > ^ij{t) • 
Q(t), let Sij{t) = Sij{t - 1), else Sij{t) = Xij{t), the global 
solution is the union of Sij (t) from all super-subSquares. 

The pseudo-codes are given in Algorithm L Since every 
node knows the locality from which it will collect information, 
it then participates the corresponding local computation, and 
at last it sends (if it is a coordinator) or receives (if not) 
the results. It may require multihop propagation to collect 
and broadcast the needed information inside super-subSquares. 
The details are omitted here. Please refer to [27 1 for de- 
tailed implementation. The coordinator computes a new local 
scheduling link set by enumeration in constant time 21^^^!, 
because the size of interference-free links for a sub-square (z, j) 
is bounded by a constant (we claim it in Lemma 1). 

Algorithm 1 Distributed Scheduling by node v under the 
linear power setting 

1: state = White; active = No; Coordinator = No; 

2: Calculate which cell node v resides in regarding to the 

current partition (iiT, at^bt)\ 
3: if V is the closest node to the center of super-subSquare 

then 

4: Coordinator = Yes; 
5: end if 

if Coordinator = Yes then 
Collect Q(t) and Sij{t - 1). 
Compute ^ij{t) in sub-square(z, j); 
if S^j{t - 1) • Q{t) > Xij{t) • Q{t) then 

else 



Sij{t) — P^ij{t)\ 
end if 

Broadcast RESULT (6'ij(t)) in super-subSquare(z, j); 
end if 

if state = White then 

if receive message RESULT (6'ij(t)) then 
if V G Sijit) then 

state = Red; active = Yes; 
else 

state = Black; active = No; 
end if 
end if 
end if 
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C. Theoretical analysis and proof 

We will then prove that our algorithm is correct and has a 
constant approximation ratio to the optimal capacity region. 

Given a network (F, E), supposing UZ^ is a set of disjoint 
local link sets inside for scheduling, where Zi e E and Zi D 
Zj = (j) if i for any link I e Zi, if I is activated, then 



= 



+ ^out 



where P denotes cumulative interference from all other 
activated links in the network, /-^ denotes the total interference 
from simultaneously transmitting links inside Zi and /^^^ 
denotes the total interference from transmissions outside. 

Therefore, we can do independent scheduling inside Zi 
without consideration of I^^^ from concurrent transmissions 
outside Zi, if I^^^ gets bounded by a constant, i.e.. 



iL < (1 - ^) • ^max, lout < e ■ /max, < S < 1, / G Z,. 

We let /max denote the maximum interference that the longest 
links in E can tolerant during a successful transmission, and 
^iax represent the maximum interference that an activated link 
/ can tolerant during a successful transmission. 

We will give a formal statement that /^^^ of each activated 
link / is indeed bounded by e/max in Lemma 2. However, to 
prove Lemma 2, we shall claim the following Lemma 1 in 
advance. 

Lemma 1: In the linear power setting, the number of 
interference-free links for a local link set Zi inside a square 
with a size length JR is bounded by a constant. Let OPTi 
refer to the optimal MWISL for Zi. That is. 



out 

^ I -'- m a x ' O 



\OPT^\ < 



{V2JRY 
{l-e) 



CT] 



+ 1, 



Proof: This proof is available in Appendix A. ■ 
We let \OPTi\uh denote an upper bound of the size of the 
local optimal MWISL for Zi, then we have. 

Lemma 2: In a given network (F, E) under the physical 
interference model in the linear power setting, if the Euclidean 
distance between any two disjoint local link sets is at least 
M X R, then activated links in each local link set suffer a 
bounded cumulative interference from all other activated link 
sets, i.e., for each activated link / in local link set Zi, 



hut < £ • /max, < £ < 1, 



where M is a constant, satisfying M > 



27rcr]RP 



■\OPTi\^ 



Proof: This proof is available in Appendix B. ■ 
To analyze the theoretical performance of our method, we 

first review the following definitions. 

Definition 2: (affectness [21]) The relative interference of 

link P on I is the increase caused by P in the inverse 

of the SINK at I, namely ri*{l) I^. /{Pir]\\uv\\-'^). For 

convenience, define ri{l) = 0. Let ci 



l-aU(Pir]\\uv\\-^) 

be a constant that indicates the extent to which the ambient 
noise approaches the required signal at receiver ty. The 
affectness of link I caused by a set S of links, is the sum 
of relative interference of the links in S on I, scaled by ci, or 
as{l) = ci ■ Yl n*(/). 

Definition 3: fp-signal set [21]) We define a p-signal set to 
be one where the affectness of any link is at most 1/p. Clearly, 
any ISL is a 1 -signal set. 



Lemma 3: ( 1211 ) There is a polynomial-time protocol that 
takes a p- signal set and refines into a p'- signal set, for p' > p, 
increasing the number of slots by a factor of at most 4(^)^. 

That is, if a p'- signal set of links can be scheduled simul- 
taneously in a single slot, then we can use some algorithm to 
partition a signal into at most 4(^)^ slots such that each set 
of links in every slot is a p'- signal set. This lemma implies that 
a p-signal set can be refined into at most 4(^)^ p'-signal set 
through a polynomial-time algorithm, e.g., a first-fit algorithm. 

Next in Lemma 4 we present a constant approximation 
relationship between each locally computed scheduling link 
set and its counterpart of the global optimal scheduling set. 

Lemma 4: The weight of Xij (t) has a constant approxima- 
tion ratio to the weight of the intersection set by the local link 
set Lij and the global optimal MWISL 

Proof: Though the total interference that every link / can 
tolerate at a successful transmission becomes (1 — s)l!^^^, 
we show that it does little influence on local computation of 
optimal link scheduling set for Lij. 

Normally any ISL is a 1 -signal set. Nevertheless, in order to 
keep independence of sub-squares, the locally computed ISL 
should be a p- signal set, where p is bigger than 1. That is, for 
the affectness of a normal ISL it holds that: 



as{l) = ci-Y n*(0 < ^-j- • ^ < 1, 

M ; I Z^^ I - l-a^/Pi Pi - 



(2) 



whereas the affectness of a locally computed ISL for sub- 
square must satisfy that. 



i*es, 



1 - a^/P, Pi 



(3) 

Therefore, by Lemma 3, any p- signal set can be refined 
into at most 4(^)^ p'-signal sets, where p' > p. So a normal 

^ ^ -Signal link sets at 



MWISL can be refined into 



(1-5)2 (1-5) 

most. Since the -Signal link set returned by enumeration 
is most weighted, so the locally computed link scheduling sets 
Xij (t) has a weight 

[l-ef 



-W{OPTij{t)). 



(4) 



Let S*j (t) = S* (t) n Lij denote the intersection by Lij and 
S*{t), where S*{t) is the global optimal MWISL at time slot 
t. Let W{S) denote the summed weight of links in a link set 
S. It is obvious that 

W{S:jit)) < W{OPTij{t)) < ^^^^W{Xij{t)). (5) 

Thus it proves the lemma. ■ 
Theorem 1: S{t) = [JSij{t) computed by our algorithm is 
an independent link set under the physical interference model 
in the linear power setting. The weight of S{t), i.e., W{S{t)), 
is a constant approximation of the weight of the global optimal 
MWISL with probability of at least 

Proof: The proof consists of two phases. We first prove 
that S{t) = [JSij{t) is an independent set. We next derive the 
approximation bound S{t) achieves. 
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Partition 
at time slot 
t*-l 



Fig. 2: Partitions at time slot t* — 1 and time slot 

Phase I: We rely on induction to infer that at every time 
slot S{t) is a union of disjoint activated local link sets that 
are separated by at least M cells from each other. Then we 
have that S{t) is an independent link set under the physical 
interference model by Lemma 2. The following are the details. 

For any link / G S{t), assuming / G Sij{t), the total inter- 
ference / suffers from all the other simultaneously transmitting 
links in S{t) is denoted by Is{t)- 

At time slot 0, every local activated link set Sij (0) = A^j (0) 
is kept 2M cells away from each other, so S{0) is an 
independent set by Lemma 2. 

At time slot 1, either Sij{0) or is chosen to be a 

part of S'(l). For those super-subSquares whose Sij{l) = 
Si'j'{0) n Yij{0), their distance is kept at least 2M cells 
away. For those super-subSquares whose Sij{l) = 
their distance is also kept at least 2M cells away. And the 
distance between the two kinds of link set is at least 2M — 1 
cells away. So S'(l) is an independent set. 

At some time slot t*, t* > 1, for some super- subSquares, 
Sij{t'^) consists of disjoint subsets from several different 
- 1) which fall into Fij(r), i.e., 

Sij(e) = 5(t*-i)ny.,(t*) = y - 1) ny,,(t*)}. (6) 

For instance, as illustrated in Fig. 2, link sets {h^h} 
and {/3,/4} respectively get scheduled in two differ- 
ent super-subSquares (i.e., super-subSquare(l, 2) and super- 
subSquare(2, 2)) at time slot t* — 1, the links {hjh^hjh} 
are then get scheduled in the same super-subSquare(l, 2) at 
time slot t*. 
We then let 



(7) 



for brevity. Each Si'j'{t'' — 1) is kept at least M cells away 
from each other, so is each 

Clearly, S{t*) can be divided into two separated subsets, 
one formed by some subsets of S{t* — 1), the other formed 
by newly computed Xij{t*), i.e., 

s[e)4js.,[e) (8) 



(9) 



mn, 
^ mn^pq 



Since U IJ <l>[,'^^., (t*) is a subset of S{t* - 1), it is composed 

pq i'j' 

by disjoint subsets with a mutual distance of M cells at least. 



The distance between any distinct is no less than 2M 

cells. Then we consider the distance between a disjoint subset 
of U U <^>[,'^^v(r) and a A'^n(^*). Since A'^n(^*) locates in 

pq i'j' 

sub-square(m, n), which is M cells away from the border of 
super-subSquare(m, n), the distance between a disjoint subset 
of U U ^^j'it*) and a ^mn{t*) is still no less than M cells. 

pq i'j' 

Comprehensively, S{t*) consists of disjoint subsets which are 
separated by at least M cells. 

Note that a disjoint subset of S{t^) does not equalize to a 
Sij{t^) since a Si'j'{t^ — 1) may be reserved completely in 
different super-subSquares at time slot t*. Here we denote the 
disjoint subset by V^i(r), and = U 

By Lemma 2, for each link / G where comes 

from the former part of equation (9), we have Is(t*) ^ ^max. 
Meanwhile, for each / G tpi{t^), where / G ipi{t^) comes from 

Then we have 



the later part of (9), it holds that /^(i*) < /^ax 



<4ax,v/G^(r), 



(10) 



indicating that S{t*) is an independent set. 

Next we consider situations at time slot t* + 1. Similarly, 
S{t* + 1) composes of disjoint subsets separated by no less 
than M cells. Using the same technique as at time slot t*, we 
can get that + 1) is still an independent set. 

By induction we can infer that S{t) is a union of disjoint 
activated subsets separated by M cells at least, thus it is an 
independent set. Herein we finish the first phrase of the proof. 

Phase II: We derive the approximation ratio by the pigeon- 
hole principle. Lemma 4, and Proposition L 

Let D{t) denote the link set of the removed strips at 
Partition ( iC, a^, D'^{t) represents a subset of links 
of which fall inside D{t), i.e., D*(t) = D(t)n6'*(t). Recalling 
that there are K'^ different partitions for a plane totally. If 
we tried all these partitions in a single slot, each cell(z, j) 
would appear in the "removed" strips at most 2KM times. 
Thus the weight of all D^{t) in the K'^ partitions should be 
2KMW{S%t)), i.e.. 



E 

i=0 



W{D*(t)) < 2KMW{S*{t)). 



(11) 



Since actually we only experience one partition at time 
slot t, by the pigeonhole principle the corresponding 
Partition ( iiT, a^, 6t) has a probability at least to be an 

optimal partition, the weight of whose removed links in D*(t) 
is no greater than that of other partitions. Instantly we have 
the following with probability of at least, 

k"^ — 1 

W{D:{t)) < • E ^(A*(t)) < ^W{S\t)), (12) 



i=0 



and then, 

WiuSljd)) = W{S''(t)\D%t)) > (1 



2M 



)W{S^t)). (13) 



Following Lemma 4, the following holds with a probability 

of 1/K\ 



2M 4 

{i-£^)w{s*{t))< 



{l-ey 



;W{UX,jit)). 



(14) 
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So by Proposition 1 we get that the approximation ratio for 
the optimal is 



(1-2.M )(i_^)2 



D. Time and communication complexity 

We then give a brief analysis on time complexity and 
communication complexity. 

We define the time complexity as the time units required by 
the running of Algorithm 1 inside each local area in the worst 
case |28|. It includes time units for local information collection 
and local computation time at the center node. At every time 
slot, the coordinators shall collect queue information and last 
scheduling status of all links inside the super- sub Square. After 
computation, the coordinators broadcast the scheduling results 
throughout the super-subSquares. It may require multihop 
propagation to collect and broadcast the needed information 
inside super-subSquares. We implement it as follows. 

To avoid collision, the coordinator can first compute a tree 
that determines the sequence of transmissions at each node in 
the super-subSquare based on topology information. As the 
root node it broadcasts the results to its children (or child). 
The children acquire the order of transmissions and relay this 
message as the specified order. The message is relayed until 
it reaches every node. Then the required information (or the 
scheduling results) can be collected (or broadcast) according to 
the sequence. We suppose that the the parent node can packet 
the children's messages into its own packet for collection. In 
the worst case each node has to transmit one by one at different 
mini slot, then it causes 0{nij) time complexity where nij is 
the number of nodes inside super-subSquare(z, j). Herein we 
just give a basic scheme, the time complexity may be further 
reduced with the help of other techniques. 

After all required information gathered, the local compu- 
tation complexity at the center node is 0(2l^^^l), which can 
be ignored comparing to the time for information collection. 
Thus the total time complexity is 0{nij). The communication 
complexity is the total number of basic messages transmitted 
during each scheduling in the worst case [28 1. The commu- 
nication complexity is 0(|F|), and the number of messages 
transmitted in each local area is 0{nij). 

IV. The algorithm in the Uniform Power Setting 

We now extend the framework to the uniform power set- 
ting. Instead of enumeration, we compute a MWISL of the 
candidate links for each sub- square by adopting the method 
proposed in ll24l . The reader may wonder why we do not 
use enumeration directly as we do in the linear power setting. 
The difficulty lies behind the fact that the cardinality of the 
optimal MWISL in each sub- square is no longer bounded by 
a constant in the uniform power setting, which drives M no 
longer bounded by a constant and simple enumeration time 
consuming. Except for the difference, the general structure of 
the algorithm is the same with that of the linear power setting. 
The following we will give a detailed analysis, and derive a 
logarithmic bound on throughput capacity. 

We first describe the main idea of Algorithm 1 in f24\ for 
computation of MWISL inside each sub-square, and show that 
the cardinality of the resulted link set is upper bounded by a 



constant. Given a set of links and weights associated with the 
links, the algorithm works as follows: 

Phase I: Remove every link whose associated weight is at 
most ^^^^ where K^max denotes the maximum weight among 
all links and n is the number of all given links. Let Wmin 
denote the minimum weight from the remaining links. 

Phase II: Partition the remaining links into log ^^^^^ groups 
such that the links of the i—th group Gi have weights within 
[2*t(;min, 2*+^i(;min]. For each group Gi of links, it finds an 
independent set of links among it by adopting the method in 
|32|. Totally it will get log ^™ ISLs, one for each group. 
Then it chooses the one with the maximum weight among the 
log ^^^^^ ISLs as the final solution. 

^ Wmin 

The remove operation of Phase I simply excludes those links 
with smaller weight so that we can get a smaller number of 
groups when partition is done according to weight. In that 
way, it is more likely to get a better approximation ratio to 
the optimal in practical, though the theoretical guarantee is 
not explicitly improved. We also notice that the resulted link 
scheduling set is an ISL by the method in 13211 . With the help 
of the following lemma presented in ll32ll . we will get that the 
the resulted link set is indeed with limited number of links. 

Lemma 5: ( (321 ) Consider a link / = (u^v) and a set N of 
nodes other than u whose distance from u is at most 
If link / succeeds in the presence of the interference from N, 
then 



liVl < 



^ R ^ 



The corollary below asserts a upper bound of the number 
of successfully transmitting links in a local link set with size 
JR X JR when utiHzing Algorithm 1 in ll24ll . 

Corollary 1: The cardinality of the resulted link scheduling 
set by Algorithm 1 in ll24ll is upper bounded, i.e.. 



1 



r 



Proof: In the method presented in |32|, it adds firstly the 
shortest link among all the candidates to the scheduling set. 
And the distance between any pair of nodes will be no greater 
than \/2JR. Therefore, let r be the shortest link length and 
then p = , we derive the upper bound. ■ 

We let \^ij{t)\ub denote a upper bound of the size of the 
local computed ISL for Zi by Algorithm 1 in [24J. Then we 
present the lemma below to show that in the uniform power 
setting, if any two disjoint local link sets of a given network 
are separated by a certain constant distance away, independent 
scheduling inside each local link set without consideration of 
interference outside is possible. 

Lemma 6: In a given network {V, E) under the physical 
interference model in the uniform power setting, if the Eu- 
clidean distance between any two disjoint local link sets is 
at least M x R , then activated links in each local set suffer 
negligible cumulative interference from all other activated link 
sets, i.e., for each activated link / in local link set Zi, 



<e-Irr 



< £ < 1, 
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where /max is the maximum interference that the longest Hnks 
in E can tolerant during a successful transmission, M is a 



27rr7P-|A'^,(t)|^b 
(/^-2)£lmaxi?^ 



constant, satisfying M > 

Proof: The proof is available in Appendix C. 
In light of Lemma 6, the partition strategy to enable 
distributed implementation will remain effective in the uniform 
power setting. 

Next, we reveal the relationship between the newly com- 
puted scheduling set for sub-square (i, j) and the corresponding 
S^j{t). It will help us to derive an approximation ratio of the 
global solution generated by our algorithm later. 

Lemma 7: The weight of the newly computed link schedul- 
ing set inside each sub-square(z, j), i.e., W{Xij{t)), has a 
approximation ratio of ^i^^y to the weight of the intersection 
set by the corresponding local link set Lij and the global 
optimal MWISL 

Proof: The newly computed local scheduling set for each 
sub- square should be a (1 — -signal set, so its weight is 



at least 



times the corresponding 1- signal set. Since the 



solution for MWISL with uniform power assignment proposed 
in ll24ll finds an ISL with affectness no greater than 1, it is 
obvious that 



W{S:Jt))<W{OPT,,{t))< 



4m 



(1 



;W{X,j{t)), (15) 



where /i = log \ V\ is the approximation bound of the algo- 
rithm in |24|. ■ 

Then, we will state an exact bound on the throughput 
performance of our algorithm in Theorem 2. 

Theorem 2: The union of the local computed scheduling 
link sets, i.e., S{t) = \JSij{t), is feasible under the uniform 
power setting. The weight of S{t) achieves a fraction of 
0(log \ V\) times the optimal solution. 

Proof: We first show that the resulted scheduling link 
set by our algorithm is an independent link set with uniform 
power at every time slot. Using the same technique in the proof 
of Theorem 1, we can inductively conclude that each global 
scheduling set S{t) is composed by disjoint link sets which 
are kept at least M cells away from each other. Therefore, by 
Lemma 6, we can get that for each link / G S{t), the total 
interference it receives satisfies that. 



'Sit) 



<{l-s)-ll 



£ ' /max ^ /i 



max -^max 



(16) 



It then follows that S{t) is an independent link set. 

We next prove the approximation bound on the performance 
of our algorithm. Remember that a Partition ( /C, a^, 6t) has a 
probability of to be a optimal partition. In that case, 

with the help of Lemma 6 it follows that 

W{UX.j{t)) > ^^~'^][^~^^^ W{US:,{t)) (17) 



{1-£)^{K -2M) 



W{S*{t)) (18) 



Hence, we get 



'{W{S{t)) > 



{l-ef{K -2M) 



S%t))>—. (19) 



Since fi = log \ V\ , it is clear that Algorithm 1 can achieve 
^ dogVl ^ optimal capacity region in the uniform power 

setting by Proposition 1. ■ 

The time complexity and communication complexity is the 
same as the corresponding of the linear power assignment. 

V. PERFORMANCE EVALUATION 

We do simulation experiments to verify the correctness, and 
evaluate the throughput performance of our proposed algo- 
rithms. The general setting of our experiments is as follows. 
We consider a network with 500 nodes, half of which as 
senders randomly located on a plane with size 200 x 200 units, 
the other half as receivers positioned uniformly at random 
inside disks of radius R around each of the senders. Packets 
arrive at each link independently in a Poisson process with 
the same average arrival rate A. Initially, we assign each link 
k packets where k is randomly chosen from [100,300]. The 
path loss exponent is set to be 3. 

We conduct three series of experiments to fully evaluate 
our algorithms. In the first series, we verify the correctness 
of our algorithms, that is, whether each activated link receives 
bounded interference from all other links outside the local link 
set it belongs to. The other two series of experiments focus 
on the evaluation of average throughput performance in terms 
of total backlog (the total number of unscheduled packets). 
The total backlog fluctuates slightly in a region if the arrival 
rate vector lies in the achievable capacity region of a link 
scheduling algorithm. Or it increases dramatically if the arrival 
rate vector exceeds the achievable capacity region. If the 
total backlog increases unboundedly to infinity, the network 
will become unstable. Thus we can roughly approximate 
an average supportable arrival rate through the changes of 
backlog. 

In the second series, we study how some related variables 
affect the performance of the algorithms. In the third series, 
we further study the performance efficiency of the proposed 
algorithms by comparisons with two centralized algorithms: 

1) Greedy Maximal Schedule (GMS), which is observed 
to perform throughput optimally in simulations [4]. It 
greedily adds to the scheduling set the most weighted 
link that is independent with the previously added links. 

2) Randomized Algorithms, where each time slot links are 
randomly picked, and added to current scheduling set if 
they do not conflicts with previous selected links. The 
process is repeatedly performed until no new links can 
be added. 

The reasons why we do not compare with distributed im- 
plementation of the two centralized algorithms are as follows. 
Undoubtedly, these two are distributed in nature according 
to some locality- sensitive approaches ll28l . if considering 
binary interference models. And previous works on binary 
interference models indeed implement them in a distributed 
manner, i.e., 13411 proposing a local greedy link scheduling 
algorithm and fSl making scheduling decisions by probability. 
The key point is that the effect of this category of interference 

^The results hold under any arrival process satisfying the strong law of 
large numbers, using the fluid model approach ll33l lfT2l . 
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t- Inference outside 
_ Maximum tolerable 
inference outside 



-•— Affectness inside 
"•"Totaiaffectness 
— Maximum affectness in; 



Link ID 

(a) Interference vs. links 




(b) Affectness vs. links 



Fig. 3: The interference/affectness each randomly chosen transmit- 
ting link receives with the linear power assignment, (a) shows the 
outside interference each link receives, (b) shows the inside affectness 
for each link, and the total affectness each receives from all other 
simultaneously transmitting links. 



"-Inference outside 
_ Maximum tolerable 
inference ouside 



-•-Affectness inside 
"•"Total affectness 
— Maximum affectness inside 



(a) Average interference vs. time slot (b) Average affectness vs. time slot 

Fig. 4: Average interference per transmitting link receives at different 
time slots with the linear power assignment, (a) shows the average 
outside interference per transmitting link receives, (b) shows the aver- 
age inside affectness, and the average total affectness per transmitting 
link receives from all other simultaneously transmitting links. 

is localized and binary. It is easy to keep feasibility of the 
final global scheduling set by inactivating all other links in 
the interference range of an activated link surely or w.h.p.. 
However, it is not so obvious and easy to guarantee under 
the physical interference model. It is hard to localize these 
approaches because the interference effect is global and addi- 
tive. Simply inactivating all other links in some neighborhood 
of an activated link will not help much to keep the final 
scheduling set feasible. Considering these difficulties, we focus 
on comparisons with the centralized visions. 

A. Under Linear Power Setting 

In this subsection we conduct the three series of experiments 
for the linear power assignment. 

We first plot the numerical results for the first series 
of experiments in Fig. 3, 4. For convenience, we call the 
interference/affectness each link receives from all other simul- 
taneously transmitting links outside/inside the corresponding 
super-subSquare it belongs to as outside interference/inside 
affectness. Since the maximum tolerable interference of each 
link for successful transmission is different, the corresponding 
inside and total interference each link receives is shown as 
the value of affectness for better illustration. Recall that, for 
a successful transmission, the inside and total affectness shall 
be less than 1 — e and 1 respectively. 

In this set of experiments, we set £ a small value of 0.3. For 
Fig. 3, we randomly choose a time slot and 10 transmitting 
links. As shown in Fig. 3(a), the outside interference each 
chosen link receives is much smaller than the constant bound 



elmax- From Fig. 3(b), we can further confirm that the trans- 
mitting links suffer negligible interference outside. For Fig. 4, 
we show the average outside interference per transmitting link 
receives at different time slots in Fig. 4(a), and the total and 
outside average affectness per transmitting link in Fig. 4(b). 
These results cooperate to show that the outside interference a 
transmitting receives is indeed bounded by elmax, and indicate 
the correctness of Algorithm 1. 

Next we set another set of experiments to study the perfor- 
mance of Algorithm 1 alone under different values of variables 
that may impact the average throughput performance actually. 
The variables ^ and e dominate the theoretical bound of 



M 



Algorithm 1. When e is fixed, a larger ^ implies a bigger 
fraction of the optimal capacity region, but with a smaller 
probability of to achieve it. e denotes a weighting factor 
between the interference a activated link suffers inside and 
outside the sub- square. A bigger value of e indicates a smaller 
value of M theoretically. Therefore, though under the fixed 
value of ^ a smaller e leads to a bigger fractional capacity 
region, the probability to achieve this region gets smaller 
because of the resulted bigger M. Since the running time of 
the simulation is much shorter than the time for the algorithm 
to achieve the theoretical value, the probability will impact 
the actual average throughput in our experiments. Typically, 
a small probability of may cause poorer throughput 
performance. An experiment study of the two variables are 
illustrated in Fig. 5 and Fig. 6. 

In Fig. 5 we compare the total backlog by increasing the 
feasible values of ^ where e serves as an constant. Fig. 5(a), 
5(b), 5(c) respectively shows the comparisons of total backlog 
at time slot 1000 with increasing arrival rate when e = 0.2, 
£ = 0.4, e = 0.8. Note that the value of ^ must be greater 
than 2, or the size of the sub- squares will be 0. It shall also 
be noticed that the feasible values of ^ are different when e 
varies, since £ affects the value of M. It then explains why we 
set different values for ^ for varying e. All the three figures 
show that the throughput performance generally improves as 
^ increases, which coincides with the theoretical results we 
derive. Though the relevant probability shall become smaller 
as ^ increases, there is no obvious sign shown in Fig. 5. (a), 
5. (b). We can see an obvious impact in Fig. 5(c) where ^ can 
be set bigger values. It shows the average throughput becomes 

6. 



a little worse at ^ 



7 than § 





Fig. 6: Total backlog vs. average arrival rate vs. different values of 
s at time slot 1000 in the linear power setting 

We briefly illustrate the impact of e at fixed ^ in Fig. 6. 
Though the theoretically achievable capacity region shall have 
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Average arrival rate 
(a) e = 0.2 




Average arrival rate 
(b) e = 0.4 




Average arrivafrate 
(c) e = 0.8 



Fig. 5: Total backlog vs. average arrival rate vs. different values of ^ at time slot 1000 in the linear power setting 



been greater with a smaller it seems to show contradicted 
results in Fig. 6(a). The apparent contradiction lies behind the 
probability of which becomes quite small because of a 
much bigger M caused by a smaller e. Fig. 6(b) shows a 
similar consequence caused by the crucial impact of £ both 
on the achievable capacity region and the relevant probability. 

We then focus on comparisons with the two centralized 
algorithms in Fig. 7, 8. We set £ = 0.8, ^ = 6 to conduct 
the following simulations. These simulation results indicate 
that our distributed scheduling algorithm(denoted by DS in 
figures) achieves comparable performance with the centralized 
GMS and the heuristic RA. 





Average arrival rate 

(a) Total backlog vs. average arrival rate 



(b) Zoom in of (a) 



Fig. 7: Total backlog vs. average arrival rate vs. different algorithms 
at time slot 1000 in the linear power setting 

In Fig. 7 we compare the total backlog changes of the 
three algorithms as average arrival rate increases. It shows 
that our algorithm may approximately support a maximum 
average rate no greater than 0.2, and the other two support no 
greater than 0.3. We zoom in a subgraph of the Fig. 7(a) in 
the Fig. 7(b), where the average arrival rate is in [0.16,0.18]. 
It shows that our algorithm can keep the total backlog stable 
in the region. Fig. 8 then illustrates detailed comparisons of 
achievable capacity region for the three algorithms. It shows 
that our algorithm can support a comparable traffic arrival rate 
vector with the two centralized algorithms. From the three 
subfigures we can see that our proposed algorithm can keep 
the total backlog stable at an arrival rate no greater than 0.18, 
while the counterpart of the greedy algorithm and random 
algorithm is 0.28 and 0.26. 



B. Under Uniform Power Setting 

Now we conduct the same three series of experiments for 
the uniform power assignment. 



The setting of the first series of experiments is the same as 
the linear power setting. The corresponding numerical results 
are plotted in Fig. 9, 10. From the two figures it concludes 
that the transmitting links receive much smaller interference 
outside than the theoretical bound. Therefore, it indicates that 
our algorithm is correctly applicable to the uniform power 
setting as well. 

Then we then conduct another set of experiments to study 
how ^ and £ affect the performance. 

We show the effect of ^ at different values of £ where £ = 
0.2, 0.4, 0.8 respectively in Fig. 11. Fig. 11(a), 11(b) and 11(c) 
plot the total backlog changes as increasing average arrival rate 
under different values of The trends are a little different 
from those in the linear power setting. From the three figures 
we can see that the average throughput performance gets better 
as increasing ^ under a fixed £. The decreasing probability 
of shows no obvious influence on the results. It may be 
partly caused by the algorithm for computing new scheduling s 
inside sub- squares since it selects candidate links based on 
their distance with previous selected links. Thus a larger area 
implies more links get scheduled. This improvement remits 
the influence of the decreasing probability of 



*- Inference outside 
_ Maximum tolerable 
inference outside 



-•-Affectness inside 
"0" Total affectness 
— IVIaximum affectness inside 




(a) Interference vs. links 



Time slot 

(b) Affectness vs. links 



Fig. 9: The interference/affectness each randomly chosen transmit- 
ting link receives with the uniform power assignment, (a) shows the 
outside interference each link receives, (b) shows the inside affectness 
for each link, and the total affectness each receives from all other 
simultaneously transmitting links. 

The similar trend occurs when £ increases at different values 
of ^ We give a brief illustration in Fig. 12. Both the figures 
show that a larger £ generates better performance at fixed ^ 
Since the affectness of local ISLs computed by the algorithm 
of 1 24 1 inside each sub-square may be much smaller than 1—5, 
it explains why £ has little influence on the theoretical bound 
actually. But e has much more influence on the value of M 
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(a) DS (b) GMS (c) RA 

Fig. 8: Achievable capacity region of different algorithms in the linear power setting 




Average arrival rate Average arrival rate Average arrival rate 



(a) £ = 0.2 (b) £ = 0.4 (c) e = 0.8 

Fig. 11: Total backlog vs. average arrival rate vs. different values of ^ at time slot 1000 in the uniform power setting 



•-Inference outside 
_ Maximum toierable 
inference ouside 




(a) Average interference vs. time slot (b) Average affectness vs. time slot 

Fig. 10: Average interference per transmitting link receives at 
different time slots with the uniform power assignment, (a) shows the 
average outside interference per transmitting link receives, (b) shows 
the average inside affectness, and the average total affectness per 
transmitting link receives from all other simultaneously transmitting 
links. 

and the corresponding probability. Therefore, in the uniform 
power setting, the algorithm achieves better performance with 
a larger e. 





Average arrival rate 

(a) # = 3 



Average arrival rate 



(b) 



Fig. 12: Total backlog vs. average arrival rate vs. different values of 
e at time slot 1000 in the uniform power setting 



We conduct the third set of experiments to compare with 
the two centralized algorithms on average throughput where 
we set £ = 0.9, ^ = 9- The results in Fig. 13, 14 convey 
the same kind of information as those in Fig. 7, 8. These 
graphs jointly show that Algorithm 1 still keep comparable 
throughput performance with the two centralized in terms of 
total backlog. For example, the maximum supportable average 
arrival rate of our algorithm is around 0.12, nearly a half of 
0.24 and 0.22 achieved by the two centralized algorithm. 




Average arrival rate 

(a)Total backlog vs. average arrival rate 



(b) Zoom in of (a) 



Fig. 13: Total backlog vs. average arrival rate vs. different algorithms 
at time slot 1000 in the uniform power setting 



C. Discussions 

Compared with experiment results with the linear power 
assignment, results under the uniform power setting are a 
little poorer in terms of total backlog under the same average 
arrival rate vector. It is probably because local computation 
of MWISL in the uniform power setting can just provide 
a logarithmic efficiency ratio, while it provides a constant 
efficiency ratio in the linear power setting. 
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Fig. 14: Achievable capacity region of different algorithms in the uniform power setting 



In both settings, our proposed distributed algorithms achieve 
almost half of the maximum supportable average arrival 
rate by the two centralized algorithms. These comparable 
results show that our distributed algorithms perform efficiently 
without much sacrifice of throughput performance. The lost 
throughput performance is mainly caused by these silent 
links outside the local scheduling sets. We have computed 
that the theoretical distance between these independent local 
scheduling sets is at least Md. In our simulation settings links 
are uniformly dispersed in the network, which means nearly 
1 - {2M/{K + 2M))^ fractional links are excluded to en- 
sure the theoretical distance for independent local scheduling. 
In our simulation results, it also shows that the theoretical 
distance is too large for practice since the actual interference 
outside the local scheduling sets are far smaller than maximum 
tolerable interference. 

From fluctuations of total backlog as time, we can roughly 
get convergence time of scheduling algorithms. In both set- 
tings, the total backlog begins to convergence to a stable 
region around time slot 2000 under supportable arrival rate 
vectors, e.g. Fig. 8, 14. The figures show that our localized 
algorithms have almost the same time of convergence as the 
two centralized algorithms. 

VI. Related work 

As a fundamental problem in multihop wireless networks, 
link scheduling has been extensively studied. Existing works 
in the literature on wireless scheduling consider different 
optimization measures and assume different interference mod- 
els. Though we mainly work on scheduling for throughput 
maximization under the physical interference model, we give 
a complete review on scheduling with distinct optimization 
measures and interference models. 

Scheduling for throughput maximization has been well un- 
derstood under binary interference models in literature ranging 
from capacity analysis to algorithm designs [10|, [HTl . [Il2l . 
161- 1^, P?!, |[5l|, (TSl. In their seminal work [29], Tassiulas 
and Ephremides prove that the celebrated maximum weighted 
scheduling (MWS) achieves optimal throughput region. Since 
finding a MWS is NP-hard in general interference models, 
they further propose a linear-complexity centralized algorithm 
using an approach called pick-and-compare [5 |. Subsequently, 
a variety of simpler and/or suboptimal scheduling algorithms 
are proposed to achieve full or fractional optimal throughput 
region under binary models. However, scheduling on a more 



practical physical interference model is only at the beginning. 

Distributed link scheduling subject to binary interference 
mainly approximates maximal scheduling algorithms locally, 
or uses the pick-and-compare approach. 

For the former, the works provide fractional throughput 
guarantee are presented in [6]- [9], [131. Using queue-length 
based random access approach, Lin and Rasool propose 
a constant overhead algorithm with 1/3 — e of the capacity 
region under the primary interference model and — e 

of the capacity region under the 2-hop interference model. 
Their works are further refined and studied in [7 1- [9 1 with 
throughput capacity guarantee no more than 1 /2 of the optimal 
under the general interference models. 

The pick-and-compare is extended to distributed scheduling 
firstly by Modiano et al [12\ in the primary interference 
model. However, it is difficult to be applied to the more 
general k-hop interference models and omits consideration 
of scheduling overheads. Sanghavi et al. |10| further present 
a distributed algorithm with constant overhead, which can 
gain arbitrary fraction of capacity region by utilizing local 
augmentation technique and the pick-and-compare. Their work 
is limited to the primary interference model as well. A novel 
distributed scheduling algorithm using graph partition under 
the pick-and-compare framework is proposed by Tang et al 
ifTTl under more mature but still binary interference models. 

Until recently, research works begin to cast attentions on 
link scheduling under the physical interference model, and 
achieve some progress. 

Goussevskaia et al fTT\ firstly present the NP-completeness 
proofs for the scheduling problem under the physical interfer- 
ence model. It firstly proposes algorithms with logarithmic 
approximation ratio 0{g{\E\)) under the uniform power set- 
ting without noise, where 0{g{\E\)) presents the link diversity 
among all links. Goussevskaia et al further [18| attempt 
to develop a constant approximation-ratio algorithm for the 
problem of maximum independent set of links (MISL), a 
special case of MWISL with uniform weight. However, it 
is valid only without noise as pointed in (35). Wan et al. 
finally succeed in developing a constant approximation-ratio 
algorithm for MISL with the existence of noise in [|32l . 

Under the linear power setting, by utilizing partition and 
shifting strategies, Xu et al. get a constant approximation 
algorithm for MWISL subject to physical interferences [|2Ql . 
Very recently Xu et al. |24| propose another constant ap- 
proximation algorithm for the same problem based on the 
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solution for the maximum weighted independent set of disks 
problem in |36|. They also develop a logarithmic approxi- 
mation algorithm for the uniform power setting. Chafekar et 
al. II22I provides algorithms for maximizing throughput with 
logarithmic approximation ratio in the two power settings as 
well. However, the attained bound is not relative to the original 
optimal throughput capacity, but to the optimal value by using 
slightly smaller power levels. 

Despite the main concern of this paper is on link scheduling 
for throughput maximization, we also make a review on the 
closely related problem of minimum length scheduling, which 
seeks a link schedule of minimum length that satisfies all link 
demands. A very recent work in [flU gives an overall analysis 
on link scheduling problems under the physical interference 
model from an algorithm view. It reveals that the algorithmic 
reduction from the minimum length scheduling to the through- 
put maximization scheduling is approximation-preserving. 

The NP-completeness proofs of the minimum length 
scheduling problem, and a logarithmic approximation algo- 
rithm without noise in the uniform power setting, is available 
in (171. With noise taken into consideration, Moscibroda et 
al. present a scheduling algorithm for the problem with power 
control, but without provable guarantee in [|15i . They subse- 
quently get a linear approximation bound in [[T6l . An attempt 
on a constant approximation bound with uniform power setting 
fails in [21], pointed out by Wan. et al |l3l. They f23l 
then propose a logarithmic approximation algorithm for the 
minimum length scheduling problem with power control. 

B lough et al. f25l claim that the so-called black links, 
with length exactly at the maximum transmission range of the 
sender, hinder a tighter approximation bounds for the mini- 
mum length scheduling problem. They try to get a constant 
approximation ratio by limiting scheduling of such kind of 
links. Thus they revise the algorithm GOW proposed in [fTTl 
by partitioning links according to the SNR diversity, instead of 
the length diversity. However, their revised algorithm GOW* 
can only guarantee a constant approximation ratio when the 
number of black links is bounded by a constant. They have 
also recognized the necessity of distributed implementation for 
the scheduling algorithm. Some discussions on the suitability 
of the algorithm for distributed execution are then presented. 

A distributed implementation of the centralized algorithm 
1371 is presented in [27 J. The algorithm assumes a variation of 
the physical interference model we study in this paper. Nodes 
rely on the carrier sensing mechanism to detect collisions, i.e., 
a node can successfully detect a carrier even if multiple nodes 
are concurrently transmitting. In each slot a node with unsat- 
isfied demand is selected as the controller. At least one link 
associated with the controller is scheduled in the slot. Other 
candidate links are added to the slot in iteratively. Whenever 
active nodes are added that cause previously allocated links 
to lose packets, this is detected through the carrier sensing 
mechanism and the SCREAM approach is used to notify all 
active nodes that the last scheduled links should be removed 
from the slot. Hence, this approach has to propagate the 
notification of SCREAM throughout the network for a feasible 
schedule. When nodes are distributed uniformly at random in a 
square of unit area, and that radio signal propagation obeys the 



log-distance path model with path loss exponent k > 2, this 
resulted schedule length is the same as the one computed by 
the centralized algorithm of [37 1. In comparison, our proposed 
localized algorithms do not need capability of carrier sensing 
and the theoretical results hold for arbitrary wireless networks. 

A recent work [38 1 proposes a column-generation-based 
solution for an optimal solution for the minimum length 
scheduling problem. They firstly formulate the problem into 
a linear programming problem, and then use the column- 
generation-based decomposition techniques to solve the prob- 
lem iteratively in fact. Though the column generation avoids an 
enumeration of all possible schedules, it still grows exponen- 
tially with the number of the links. The method with so high 
computing complexity, like the optimal primal-dual solution 
for binary interference models, is unrealistic for wireless 
networks in nature, though optimal. 

A closely related work ['261 analyzes the performance of 
longest queue first (LQF) algorithm under the physical model 
with the uniform power assignment, and employs a technique 
named "interference localization" to prevent the performance 
of LQF vanishing, and enable distributed implementation of 
LQF. The implementation of the distributed LQF algorithm 
actually relies on the radius of interference neighborhood. 
However, the trivial procedure for determination of interfer- 
ence neighborhood is centralized, and iterated at every link to 
ensure feasibility of the scheduling. Obviously it is impractical 
in large networks since number of links is unbounded. Further, 
the achievable throughput capacity is ^^^^ fraction of the 
original throughput where Jlqf is efficiency ratio of LQF 
under interference localization constraints. Despite the value 
of JLQF is bounded away from zero, but still uncertain. 
In contrast, the efficiency ratio of our algorithm under the 
uniform power setting is logarithmic. Meanwhile, our algo- 
rithm is much simpler and more practical. It implements 
distributed scheduling by compulsory separation of local link 
sets where the separation distance is predefined. We do not 
need globally collected information in the scheduling process, 
but just some neighborhood information in constant distance. 
Moreover, we consider links transmits at heterogeneous power 
and aim at distributed throughput-optimum scheduling, but not 
improvement of a specified scheduling algorithm LQF. 

Another closely related work in [39] proposes a CSMA- 
type distributed link scheduling approach, where a link that 
is not active at time slot t — 1 randomly waits for a backoff 
time to be active if it does make links in global scheduling 
S{t — 1) violate the SINR threshold. In order to compute 
a feasible schedule, this random access based approach has 
to consume much time on the backoff process. Moreover, 
it induces high communication overhead since whenever a 
candidate link is to be activated, the corresponding notification 
should reach all previously added links for verification of the 
SINR requirement. Conversely, our algorithm computes local 
results parallelly with local information within constant hops. 

VII. Discussions 

Note that our proposed algorithms have assumed a static 
wireless network where every node knows its own location and 



14 



locations of neighbors within a distance of R. It is necessary 
in the partition process and in the Hnear power assignment. To 
get all these required location information, we only need to 
perform some localization techniques with high accuracy once 
for a static wireless network, e.g., storing all these information 
at each node when deploying networks. Thus we do not 
necessarily need an extra location device, e.g., GPS device, 
embedded in every node. As to the impact of some inaccurate 
locations on the performance of our algorithms, we argue that 

1) It has little influence on the correctness of the al- 
gorithms. The distance we set for separation of each 
local area is safe enough to ensure local independent 
scheduling, since we have considered the worst case 
(with a maximum power and a highest node density 
that a feasible schedule can tolerate) when computing 
the distance. The distance can be set even larger if 
considering inaccurate locations. An inaccurate location 
may influence the decision that whether a node shall 
participate in a local scheduling, affecting the local 
scheduling set but not the correctness of the algorithms. 

2) Inaccurate locations will impact local scheduling sets, 
and thus the throughput performance. An inaccurate 
location may lead to an inaccurate link length. Since 
the link length is bounded by a minimum length and 
a maximum length both previously known, the error 
link length shall also lies in this region. The assigned 
powers in the linear power setting remain in the original 
region, though some links have a power bigger than 
they should have and some links have a smaller power. 
Consequently, the size and feasibility of local scheduling 
sets will be partly affected. To ensure the feasibility 
of local scheduling set, we can deliberately assume 
a smaller maximum tolerable inside interference for 
each link. The size and then weight of the resulted 
local scheduling sets might be smaller than the optimal 
local scheduling sets. From simulation results, the actual 
inside interference suffered by each activated link is 
much smaller than the theoretical value. Thus we infer 
that the influence won't be too serious if just some 
inaccurate locations occurs. 

We have to predefine e, K and M to initiate the algorithm 
for the partition and shifting process. But we do not need 
other extra information to determine the values, except the 
topology information. We just need to set e among (0, 1), and 
a greater value leads to better performance. We shall set M 
as some constant no less than the lower bound and K/M 
greater than 2. The exact values of K and M also depend 
on the number of squares we want to partition the network. 
Combining our simulations, it is better to set K/M around 6 
for better throughput performance, but not necessarily. As to 
the value of longest and shortest link length, these values are 
also easy to get as long as the topology is determined. We 
shall notice that once the longest link length R is determined 
as a constant, /max is also determined as a constant, even for 
the case when /max is close to 0. In such a case, the optimal 
capacity region of the network itself is also constrained by the 
small value of /max- Thus we do not have to take extra care of 



long links. The passive impact of long links on capacity region 
instead suggests that we shall avoid long links in network 
deployment if possible. 

Despite we only consider single hop link scheduling at 
each link, our work is extendible to the joint routing and 
scheduling problem using the back-pressure mechanism [33 1. 
The difference is that the weight of each link becomes the 
difference of unscheduled packets destined to the same node 
at two neighboring nodes. 

Throughout the paper we have mainly focused on link 
scheduling for optimal throughput. Another important opti- 
mization objective for scheduling is delay. It is intuitive that 
MWM-type scheduling which tries to minimize total queue 
length may suffer serious delay and long convergence time, 
detailed analysis seen in |40| [41] and their references. In gen- 
eral MWM-type scheduling, factors such as time complexity 
of scheduling algorithms, duration of data transmission and 
infrequent schedule updates impact delay. Some theoretical 
analysis and engineering implications is available ll42ll . Many 
techniques such as "shadow queues" |40|, or cross-layer 
methods [43 ], are proposed to improve delay performance. We 
leave it as future works to study efficient localized scheduling 
algorithms for different delay and throughput demands. 

VIII. Conclusion 

We tackle the problem of throughput-optimum localized link 
scheduling in multihop wireless networks subject to physical 
interference constraints. We successfully address the primary 
challenge of global interference constraint in the physical 
interference model that hinders the development of localized 
scheduling algorithms. We then utilize the graph partition and 
shifting techniques to design a localized scheduling method 
with a constant approximation ratio to the optimal solution 
under the linear power setting. We further extend the solutions 
to the uniform power setting, and develop a localized algo- 
rithm with a logarithmic approximation ratio. Our extensive 
simulation results show that under both power settings the 
achieved capacity regions of our algorithms are nearly half of 
these achieved by counterpart centralized algorithms, which 
further demonstrates performance efficiency of our algorithms. 
We believe that our work can find applications in some time- 
slotted wireless networks, e.g., time-slotted wireless sensor 
networks or wireless mesh networks. Knowing these factors 
that determine the throughput of a network also helps to better 
deploy a multihop wireless network and enhance the overall 
throughput performance. 

For future work, there remain many interesting open ques- 
tions. The first is to further improve the efficiency of localized 
link scheduling algorithms. The second is to improve the 
approximation ratio of a method in uniform power setting. 
Up to now there have been only logarithmic approximation 
algorithms with the uniform power assignment. In addition, 
our results as well as many other results, assume the whole 
topology and geometry location of the network known by each 
node. This information is not always available and accurate 
in some cases. Some simple distributed methods without 
using explicit global and location information may be more 
interesting and practical. 
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Appendix A 
Proof of Lemma 1 

Proof: Suppose T is the shortest hnk in OPTi, let I 
denote any other Hnk different from /* in OPTi. Assuming 
/* = / = {u,v), we have \\uv*\\ < V2JR. The relative 

interference /* received from OPTi is shown below. 
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Appendix B 
Proof of Lemma 2 

Proof: For any two links /* and / which do not belong to 
the same local link set Zi, assuming /* = (i^*, ^'*), / = {u, v), 
it always holds that \\uv''\\ > M x R . 

The total interference T suffered form all concurrent trans- 
mitting links located in other sub- squares is 
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since is a constant greater than 2 typically. 
To assure lout < ^ • ^max, we derive that 
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where /max = """^^^ - ^ is a constant for a given network, 
\OPTi\ub is a constant upper bound of \OPTi\. ■ 



Appendix C 
Proof of Lemma 6 

Proof: For any two links /* and / which do not belong to 
the same local link set Zi, assuming /* = (i^*, t'*), I = [u, v), 
it always holds that > M x R . The total interference 

/* suffered form all concurrent transmitting links located in 
other sub- squares is 

Ct = E pxv>^\\uv*\r 

ie{uZi\z'},i*ez' 
< Pxr^x ^^^^ X \^^j{t)U X {xMR)--dx 
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since is a constant greater than 2 typically. 
To assure I^^^ < s-Ima^, we derive M > 
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where /max = ^ is a constant, and is a 

constant upper bound of ■ 



